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Intelligent,  Robust  Control  of  Deteriorated  Turbofan  Engines  via 
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Jonathan  S.  Litt 
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A  method  for  accommodating  engine  deterioration  via  a  scheduled  Linear  Parameter  Varying  Quadratic 
Lyapunov  Function  (LPVQLF)-Based  controller  is  presented.  The  LPVQLF  design  methodology  provides  a 
means  for  developing  unconditionally  stable,  robust  control  of  Linear  Parameter  Varying  (LPV)  systems.  The 
controller  is  scheduled  on  the  Engine  Deterioration  Index,  a  function  of  estimated  parameters  that  relate  to 
engine  health,  and  is  computed  using  a  multilayer  feedforward  neural  network.  Acceptable  thrust  response 
and  tight  control  of  exhaust  gas  temperature  (EGT)  is  accomplished  by  adjusting  the  performance  weights  on 
these  parameters  for  different  levels  of  engine  degradation.  Nonlinear  simulations  demonstrate  that  the 
controller  achieves  specified  performance  objectives  while  being  robust  to  engine  deterioration  as  well  as 
engine-to-engine  variations. 


Nomenclature 


A,  B,  C,  D 

CLP 

5 

A 

G 

FN,  SM,  T5 
K 

L,M 

P 

W 

x,  y,  p,  u 
v 
w 
z 


State  Space  Model  System  Matrices 
Subscript  denoting  closed-loop  system 
Deviation  from  steady  state  operating  condition 
Block  diagram  representation  of  model  uncertainty 
Block  diagram  representation  of  plant 

Engine  net  thrust,  stall  margin,  and  exhaust  gas  temperature  respectively 

Subscript  denoting  controller,  block  diagram  representation  of  controller 

State  Space  health  parameter  augmentation  matrices 

Lyapunov  Function  matrix  for  unforced  system,  scheduling  parameter  space 

Lyapunov  Function  matrix  for  forced  system 

State  variable,  system  output,  health  parameter,  and  system  input. 

Controller  inputs 
Exogenous  system  inputs 

Unmeasured  system  outputs  related  to  control  system  performance 


I.  Introduction 

Turbofan  engine  performance  varies  from  engine  to  engine  due  to  manufacturing  tolerances,  aging,  and  deterioration 
caused  by  use.1  Generally  the  control  system  developed  for  the  engine  is  robust  enough  to  keep  it  operating  within 
acceptable  boundaries  for  several  thousand  flight  cycles,  even  though  the  degradation  will  eventually  require  the 
engine  to  be  overhauled  as  limits  are  reached.  These  limits  include  operability  constraints  such  as  maximum 
temperatures,  minimum  stall  margins,  and  performance  constraints  such  as  the  FAA’s  rise  time  requirement  for 
thrust  in  commercial  engines.  Exhaust  gas  temperature  (EGT)  increases  as  the  engine  ages  and  is  used  as  an 
indicator  of  the  engine’s  internal  condition,  i.e.,  unmeasurable  temperatures  that  can  lead  to  thermomechanical 
fatigue  (TMF)  of  engine  components.  EGT  is  usually  at  a  maximum  during  the  climb  phase  of  flight,  and  if  it 
reaches  a  specified  threshold  it  signifies  that  the  engine  needs  to  be  overhauled.  Stall  margin,  or  the  distance  an 
engine  component  is  operating  from  the  stall  line  on  its  performance  map,  tends  to  decrease  with  age,  increasing  the 
risk  of  an  engine  stall  during  transient  operation. 


NASA/TM— 2004-213375 


1 


Turbofan  engines  typically  control  Engine  Pressure  Ratio  (EPR)  or  fan  speed  to  generate  the  desired  thrust,  since 
thrust  can  not  be  measured  directly  during  flight.  As  the  engine  is  used,  manifestations  of  wear  such  as  turbine  blade 
erosion,  larger  clearances,  etc.,  degrade  the  engine’s  performance.  In  order  to  achieve  the  same  level  of  thrust  as  in  a 
new  engine,  a  deteriorated  engine  must  run  hotter  and/or  faster.  This  shift  from  nominal  operation  increases  with 
use,  and  eventually  reaches  the  point  where  the  engine’s  performance  can  not  be  maintained  without  compromising 
the  life  of  its  components.  While  regulated  variables  are  maintained  at  their  set  points  regardless  of  engine 
degradation,  the  non-regulated  parameters  shift  from  their  nominal  values  with  deterioration.2  Additionally,  in  the 
degraded  engine,  the  actual  thrust  output,  which  is  indirectly  controlled  through  the  regulation  of  other  variables, 
may  be  shifted  from  the  expected  value. 

In  this  work  the  objective  is  to  maintain  thrust,  even  in  the  presence  of  degradation,  which  precludes  the  use  of  EPR 
or  fan  speed  as  the  controlled  variable.  Here  thrust  is  controlled  directly,  given  the  assumption  that  a  reliable 
estimate  is  available.  Additionally  the  controller  minimizes  exhaust  gas  temperature,  in  terms  of  the  steady  state 
value  and  in  excursions  during  transients,  while  maintaining  acceptable  component  (estimated)  stall  margins  to 
overcome  the  operability  problems  associated  with  traditional  controllers  applied  to  aging  engines.  To  meet  these 
goals,  a  robust  Linear  Parameter  Varying  (LPV)  controller3  is  designed  to  accommodate  a  range  of  degradation 
using  sets  of  gains  scheduled  on  the  engine  degradation  index  (EDI).  The  Linear  Parameter  Varing  Quadratic 
Lyapunov  Function  (LPVQLF)  design  methodology  provides  a  means  for  developing  unconditionally  stable,  robust 
gain  scheduled  control  of  nonlinear  plants  with  appropriate  Jacobian  linearizations  across  an  operating  envelope  i.e., 
Linear  Parameter  Varying  (LPV)  systems.  The  LPV  controller  is  applied  to  a  nonlinear  simulation  of  a  large 
commercial  turbofan  engine.  This  paper  is  organized  as  follows:  section  II  describes  the  characteristics  of  engine 
degradation  and  its  accomodation  using  LPV  control,  section  III  presents  an  overview  of  the  LPV  controller 
development  process  via  Quadratic  Lyapunov  Functions,  section  IV  provides  the  pertinant  details  concerning  the 
design  of  LPV  thrust/EGT  controllers  for  a  large  turbofan  engine,  section  V  presents  the  results  of  nonlinear 
simulation  testing  of  these  controllers  and  section  VI  presents  conclusions  and  future  work. 

II.  Engine  Degradation  and  LPV  Control 

To  develop  the  LPV  controller,  models  of  the  engine  ranging  from  nominal  to  fully  degraded  need  to  be  generated. 
Off-nominal  values  of  specific  internal  engine  parameters  representing  component  efficiencies  and  flow  capacities 
are  often  used  to  account  for  these  performance  variations.  These  adjustment  parameters  are  called 
health  parameters  4  because  they  indicate  the  level  of  engine  deterioration.  The  equations  describing  the  degraded 
engine’s  behavior  are  given  by 


X(t)  =  f(x(t),u(t),p)  (2.1) 

y(t)  =  g  (x(t),u(t),  p) 

where  p  represents  the  vector  of  health  parameters.  When  obtaining  a  standard  linear  point  model  of  an  engine,  the 
health  parameters  are  treated  like  inputs. 


8x(t)  =  A8x(t)  +  B8u(t)  +  L8p  ^  2) 

8y(t)  =  C8x(t)  +  D8u(t)  +  M  8p 

Depending  upon  how  the  health  parameters  manifest  themselves,  the  system  dynamics  may  or  may  not  change  with 
degradation.  Equation  2.2  clearly  illustrates  that  steady  state  is  only  obtained  when  the  x(t)  and  u(t)  vectors  shift  to 
compensate  for  p,  and  the  output  equation  shows  how  nonzero  values  of  p  can  produce  additional  steady  state  shifts 
in  the  output  variables.  These  equations  also  imply  that  degradation  causes  shifts  in  the  engine’s  trim  values,  and  it 
is  these  shifts  that  can  result  in  unacceptable  operation.  Figure  2.1  presents  the  effect  of  engine  deterioration  on 
exhaust  gas  temperature  (EGT)  open  loop  response  (deviation  from  nominal)  to  a  2%  change  in  fuel  flow.  A  shift  in 
trim  is  clearly  evident,  with  a  less  pronounced  change  in  EGT  dynamic  response. 
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Figure  2.1. — Open  loop  EGT  response  (deviation  from  nominal)  to  a  2%  increase  in 
fuel  flow  for  several  levels  of  deterioration. 

In  general,  the  health  parameters  vary  slowly  enough  with  time  that  they  are  treated  as  constants  in  equation  2.2. 
Trending  of  the  health  parameters  can  be  achieved  by  an  estimation  algorithm  known  as  a  tracking  filter.  Tracking 
filters  have  been  successfully  used  to  estimate  health  parameters  from  measured  variables,5  allowing  more  accurate 
tracking  of  unmeasured  engine  variables  such  as  thrust  and  stall  margin.  In  this  work  it  is  assumed  that  reliable 
estimates  of  the  health  parameters,  thrust  and  stall  margin  are  available. 

As  mentioned  previously,  the  LPV  controller  is  scheduled  on  EDI.  Here  the  health  parameters  are  each  assumed  to 
change  as  functions  of  EDI,  with  some  variation  about  a  nominal  degradation  profile.  For  the  present  study,  a 
multilayer  perceptron  neural  network  is  used  to  map  estimated  engine  health  parameter  values  (assumed  to  be 
provided  by  a  tracking  filter)  back  to  the  scalar  EDI  to  accommodate  uncertainties  associated  with  the  health 
parameters.  Desired  response  of  engine  thrust  and  exhaust  gas  temperature  (EGT)  is  accomplished  by  adjusting  the 
performance  weights  on  these  parameters  for  different  levels  of  engine  degradation.  Model  uncertainty  (i.e.,  engine- 
to-engine  variation)  is  also  incorporated  into  the  design.  Figure  2.2  presents  an  overview  of  the  deterioration-tolerant 
LPV  controller. 
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Figure  2.2. — Overview  of  turbofan  engine  deterioration-tolerant  LPV  controller. 

III.  Linear  Parameter  Varying  Controllers  Based  on  Quadratic  Lyapunov  Functions 

Unconditional  stability  of  parameter  dependant  controllers  is  provided  via  development  and  solution  of  a  suitable 
Lyapunov  equation,  which  incorporates  the  dynamic  characteristics  of  the  closed  loop  system.  The  Linear  Parameter 
Varying  (LPV)  system  considered  here  is  such  that  all  time  dependence  in  the  system  is  correlated  to  the  parameter 
dependence  on  time,  with  each  set  of  state  space  matrices  used  for  controller  design  a  function  of  the  EDI,  p(t). 

Consider  first  an  unforced,  time-invariant  nonlinear  system  of  the  form 

x(t)  =  f(x(t))  (3.1) 

with  origin  x(0).  The  above  is  Lyapunov  stable  if  a  positive  definite  function,  V(x),  exists  such  that 

V(x)  =  ^Wf(x)  (3-2) 

dx 

is  negative  semi-definite.  If  equation  3.2  is  negative  definite,  the  origin  is  asymptotically  stable.  For  linear  time 
invariant  (LTI)  systems 

x(t)  =  Ax(t)  (3.3) 

and  a  candidate  Lyapunov  function  takes  the  form 


V(x)  =  xTPx 


(3.4) 


with 


V(x)  =  xT  (atP  +  Pa)x 


(3.5) 
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If  a  matrix  P  can  be  found  such  that  (AT  P  +  PA)  is  negative  (semi)  definite  then  the  system  is  asymptotically 
(Lyapunov)  stable.6  For  forced  LTI  systems 

x(t)  =  Ax(t)  +  Bu(t)  p  6) 

y(t)  =  Cx(t) 

the  Bounded  Real  Lemma6  provides  a  Lyapunov-like  stability  criterion  i.e.,  the  system  is  internally  exponentially 
stable  and 


|c(si-a)-'b|L 


if  and  only  if  the  matrix  X  =  XT  >  0  provides  a  solution  to 

AtX + XA+  XBEF  X + CTC  <  0  (3.8) 

Consider  a  nonlinear  system  approximated  by  an  LPV  system,  where  p(t)  e  P,  P  being  the  space  of  all  possible 
parameter  trajectories, 


x(t)|  A(p(t))  Bw(p(t))  Bu(p(t))  x(t) 

z(t)  =  Cz(p(t))  Dzw(p(t))  Dz(p(t))  w(t) 

y(0  c(p(t))  Dw(P(t))  o  U(t) 


with  the  feedback  controller  defined  as 


jAK(p(t))  Bk 
CK(p(t))  Dk 


(p(0)TxK(t) 
(P(t»  L  y(t) 


The  corresponding  closed-loop  LPV  system  is  represented  by 


^ clp (0  _  ACLP(p(t))  BCLP(p(t))  xCLP(t) 
z(t)  J  LCclp(p(0)  DCLP(p(t))J[  w(t) 


(3.11) 


where 


xClp  =  [xt  x^f 


Aclp(p)  - 


"  A(p )  +  B(p  )DK  (p  )C  B(p)C 
BK(p)C(p)  Ak( 


)CK(p)' 

x(P) 


®CLp(P)  — 


Bw  (p  )  +  B(p)DK  (p  )Dw  (p  )" 
BK(p)Dw(p) 


(3.14) 


Cclp(p)  =  [Cz(p)  +  Dz(p)DK(p)C(p)  Dz(p)CK(p)] 
DCLP(p)  =  [Dzw(p)  +  Dz(p)DK(p)Dw(p)] 


with  the  time  dependence  of  p  implied.  If  the  matrix  inequality 
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AjLp(p)  W(P)  +  W(p)ACLP(p)  W(p)BCLP(p) 

BclPw(p)  -yi 

^CLp(P)  BCLP(p) 

and  its  Lyapunov  solution  matrix 

W(p)  >  0 

for  all  permissible  trajectories  of  the  scheduling  parameter  p  e  P,  then  a  family  of  controllers  dependant  on  p 
(equation  3.10  )  will  provide  globally  stable,  induced  h2  norm  performance  in  terms  of  y7,8  i.e., 

II  4<y2H  <3-19) 

Equation  3.17  is  nonlinear  and  thus,  may  not  be  convex.  However  with  W  independent  of  p,  equation  3.17  provides 
a  Linear  Matrix  Inequality  (LMI)  allowing  for  numerical  solution  via  state-of-the-art  semidefinite  programming 
(SDP)  techniques.7,8  For  an  LTI  system,  this  reduces  to  determining  the  y-bounded  H„  norm.  The  controller 
formulas  (equation  3.10)  are  explicit  functions  of  feasible  solutions  to  additional  LMI-based  constraints.7,8  For  the 
present  study,  the  parameter  space  P  (i.e.,  the  range  of  EDI  under  consideration)  is  discretized  with  the  solution  of 
equation  3.17  constrained  to  be  a  single  Lyapunov  matrix  solution  for  all  of  the  corresponding  plants. 

IV.  Controller  Development 

LPVQLF  control  seeks  to  provide  a  globally  stable  closed-loop  system  with  optimal  performance  described  by 
Equation  3.19.  This  is  achieved  by  accommodating  uncertainty  in  the  system  via  1)  the  parameter  dependence 
(deterioration  dependence)  of  the  state  space  system  described  by  equation  3.9,  and  2)  the  ability  to  incorporate 
model  uncertainty  at  each  deterioration  point  (i.e.,  engine  to  engine  variation).  The  interconnection  block  diagram 
used  for  LPV  control  system  design  is  presented  in  Figure  4.1.  The  nominal  plant  (i.e.,  the  plant  used  for  controller 
design)  is  shown  as  G  and  is  the  state-space  model  of  a  large  turbofan  engine.  Gp  is  the  perturbed  plant,  which 
includes  the  effects  of  uncertainty  due  to  model  inaccuracies  (e.g.,  engine-to-engine  variation),  and  Ga  is  the 
augmented  plant  incorporating  performance  weights  for  controller  design.  Plant  inputs  are  specified  by  w 
(exogenous)  and  u  (control).  Plant  outputs  are  specified  by  z  (unmeasured  performance  measures)  and  y 
(measurements).  For  the  large  commercial  turbofan  engine  simulation  used  in  this  work,  the  available  actuators  are 
fuel  flow  (WF),  variable  bleed  valve  position  (VBV),  and  variable  stator  vane  (VSV)  position.  An  exhaust  gas 
temperature  measurement,  as  well  as  accurate  thrust  and  stall  margin  estimates  are  assumed  to  be  available. 

Optimal  performance  is  achieved  by  specifying  a  set  of  performance  weights  included  in  the  augmented  plant,  Ga,  of 
the  system  along  with  the  nominal  plant  G.  By  satisfying  equation  3.19  for  a  specified  value  of  y  of  the  augmented 
plant,  the  induced  L2  norm  of  the  disturbance-to-performance  measure  transfer  matrix  is  maintained  below  y  for  all 
plants  specified  by  p  e  P.  Robustness  to  plant  uncertainty  (i.e.,  in-band  model  error  or  neglected  dynamics)  and 
minimization  of  input  control  energy  are  also  addressed  by  including  their  effects  or  allowed  limits  via  suitably 
designed  frequency-dependant  weights  (e.g.,  Wfn  and  Wy).  The  goals  of  the  controller  are  to  provide  acceptable 
thrust  response  while  minimizing  the  net  change  in  exhaust  gas  temperature  in  the  presence  of  engine  deterioration 
and  model  uncertainty.  This  is  analogous  to  the  so-called  mixed  sensitivity  minimization  problem  solved  via  H„- 
optimal  control.9 

Shown  in  Figure  4.1  as  an  exogenous  input  to  the  plant,  G,  are  norm-bounded  perturbations  in  the  plant  model  i.e., 
uncertainty.  Since  a  detailed  uncertainty  representation  may  add  significantly  to  the  complexity  of  the  controller 
(i.e.,  controller  order),  an  unstructured  diagonal  input  uncertainty  model  is  used  with  the  upper  bound  set  as  the 
frequency  dependent  weight  Wi(s).  A  simplified  transfer  function  is  assigned  to  envelop  the  frequency-dependant 
error  between  the  “actual”  and  reduced  order  plants.  Wi(s)  essentially  represents  the  normalized  difference  between 
the  actual  (or  perturbed)  and  the  nominal  (reduced  order)  plant  i.e., 


Cclp(p) 

Dclp(P) 

— yi 


<o 


(3.17) 


(3.18) 


NAS  A/TM— 2004-2 13375 


6 


To  reduce  the  order  of  the  controller,  the  frequency-dependant  weights  were  set  to  constants,  effectively 
representing  300%  modeling  error  across  the  entire  frequency  range  of  interest.  Plant  interconnection  and  controller 
design  was  performed  using  the  MATHWORKS  Mu  Analysis  and  Synthesis  and  LMI  Toolboxes  as  part  of  the 
MATLAB  software  package.10’11 


v 

Figure  4. 1 . — Interconnection  diagram  for 
LPVQLF-based  turbofan  engine  thrust/EGT  controller. 


The  deterioration  index,  itself  being  a  function  of  engine  health  parameters  such  as  component  efficiencies  and  flow 
capacities,  provides  the  controller  scheduling  parameter  p.  The  uncertainty  associated  with  the  deterioration  index  is 
primarily  due  to  the  estimation  of  the  health  parameters  i.e.,  each  deterioration  index  is  associated  with  a  specific  set 
of  health  parameters.  However,  in  service  the  likelihood  of  a  specific  set  of  health  parameters  corresponding  to  one 
of  the  parameter-dependant  plants  used  in  the  controller  design  may  be  extremely  low.  Thus,  a  mapping  from  health 
parameters  to  deterioration  index,  capable  of  accommodating  slight  variations  in  health  parameters  from  the  values 
specifically  used  in  the  controller  design,  is  required  for  calculation  of  the  scheduling  parameter  p.  As  an  initial 
approach,  the  simple  feedforward  neural  network  shown  in  Figure  4.2  was  chosen  to  map  the  estimated  (uncertain) 
health  parameter  set  to  the  scheduling  parameter,  EDI. 
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The  set  of  sixteen  health  parameters,  known  to  vary  with  engine  deterioration  and  engine  to  engine  variation,  is  fed 
through  input  layer  weights  to  hidden  layer  neurons,  which  provide  an  output  based  on  the  degree  of  activation  of 
Sigmoidal  transfer  functions.  The  neuron’s  weighted  outputs  are  fed  to  the  single  output  layer  neuron,  whose  linear 
activation  function  determines  the  network’s  final  output  i.e.,  the  EDI.  A  training  set  consists  of  an  expected  set  of 
health  parameters  corresponding  to  a  specified  level  of  deterioration,  which  assumes  an  overall  correspondence 
among  the  health  parameters  with  respect  to  engine  deterioration.  The  network  was  trained  on  health  parameter  sets 
corresponding  to  0%,  10%,  20%,  30%,  40%,  50%,  60%,  70%,  and  80%  levels  of  deterioration  using  the 
MATHWORKS  Neural  Network  Toolbox  as  part  of  the  MATLAB  software  package.12  Refinement  of  the  training 
sets  (e.g.,  additional  training  sets  and  adding  slight  amounts  of  noise)  may  be  accomplished  in  the  future  to  enhance 
the  functionality  of  the  network. 


Figure  4.2. — Feedforward  neural  network  used  to  map  engine  health  parameters  to 
the  controller  scheduling  parameter,  EDI. 


V.  Nonlinear  Simulation  Testing 

The  following  presents  the  results  of  testing  the  LPV  controller  on  a  nonlinear  simulation  of  a  large  turbofan  engine. 
Nine  interconnected  plants,  corresponding  to  0%,  10%,  20%,  30%,  40%,  50%,  60%,  70%,  and  80%  levels  of 
deterioration,  were  provided  to  an  LMI-based  convex  optimization  algorithm,6  which  searches  for  a  feasible  solution 
to  equation  3.17  through  3.19  -  the  output  being  the  Lyapunov  function  matrix,  equation  3.17.  Additional  LMI 
contraints  incorporated  into  the  optimization  provide  the  controller  state-space  matrices  AK,  BK,  CK,  and  DK.8 

Figures  5.1,  5.2,  and  5.3  present  a  comparison  of  a  thrust-only  controller,  (EGT  varying  uncontrolled)  and  a 
controller  incorporating  EGT  control  as  well  as  thrust  control.  The  nonlinear  engine  model  is  40%  deteriorated 
(deterioration  index  of  0.4)  at  standard  day  cruise  conditions,  with  the  expected  set  of  health  parameters 
corresponding  to  this  level  of  deterioration.  Thrust  is  controlled  primarily  via  the  fuel  flow,  WF,  while  tight  control 
of  EGT  requires  manipulation  of  VBV  and  VSV  as  well.  For  the  non-EGT  controlled  case  (fig.  5.1),  a  3%  increase 
in  thrust  demand  results  in  an  11.5  degree  increase  in  EGT  (fig.  5.1b).  VSV  and  VBV  are  unactuated  for  this 
situation  (fig.  5.1c).  In  contrast,  for  the  same  thrust  demand  input,  the  EGT-controlled  engine  (fig.  5.2)  attains  a 
increase  in  EGT  of  7.6  degrees  -  a  net  reduction  of  3.9  degrees  in  EGT  relative  to  the  uncontrolled  case.  Figure  5.2c 
presents  the  corresponding  changes  in  VSV  and  VBV  for  the  transient,  illustrating  the  effectiveness  of  these 
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actuators  in  EGT  control.  The  positions  of  these  actuators  are  typically  scheduled  on  engine  operating  conditions, 
thus  the  deviations  presented  would  be  added  to  the  positions  dictated  by  the  schedule.  Figure  5.3  presents  the 
corresponding  component  stall  margin  deviations  (SMD)  for  the  non-EGT  controlled  and  EGT  controlled  cases.  As 
shown  in  the  figure,  the  SMD  show  no  significant  change  -  all  values  are  well  above  the  stall  limits  for  each 
component. 


0  10  20  30  00  SO  60  70  80 


Time  («e) 

Figure  5. 1 . — Thrust-only  control  of  40%  deteriorated  engine  (EDI  of  0.4). 
a)  Normalized  thrust  response,  b)  EGT  deviation  from  steady  state, 
c)  VSV  and  VBV  deviations  from  steady  state. 


0  10  20  30  40  50  00  70  00 
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Figure  5.2. — Thrust  and  EGT  control  of  40%  deteriorated  engine  (EDI  of  0.4). 
a)  Normalized  thrust  response,  b)  EGT  deviation  from  steady  state, 
c)  VSV  and  VBV  deviations  from  steady  state. 
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Figure  5.3. — Stall  Margin  Deviation  (SMD)  response  of  thrust  and  EGT 
control  of  40%  deteriorated  engine  (EDI  of  0.4). 
a)Thrust-only  control,  b)  Thrust  and  EGT  control. 


As  mentioned  previously,  the  parameter-dependant  controllers  are  designed  to  accommodate  plant  uncertainty  due 
to  engine  deterioration  (i.e,  controller  scheduled  on  EDI),  with  considerable  uncertainty  incorporated  into  the  design 
due  to  engine-to-engine  variation.  Other  control  methods  e.g.,  Hx  control,9  while  accommodating  plant  uncertainty 
at  a  specific  operating  point  in  a  similar  fashion,  are  intended  for  LTI  systems  and  fail  to  adequately  accommodate 
modeling  errors  due  to  significant  plant  nonlinearities  across  an  operational  envelope.  Consider  the  LPV-controlled 
engine  response  shown  in  Figure  5.4.  For  the  nominal  engine,  the  corresponding  LPV  controller,  robust  to  engine- 
to-engine  variations,  achieves  the  thrust/EGT  response  objectives. 


t 
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Figure  5.4. — Nominal  (undeteriorated)  engine  response  with  LPV 
controller  scheduling  parameter  (EDI)  equal  to  0.0. 
a)  Thrust  response,  b)  EGT  deviation  response,  c)  VSV  and  VBV  response. 


Figure  5.5  presents  a  70%  degraded  engine  controlled  by  the  nominal  engine  controller  (thrust  and  EGT  controlled). 
Although  the  thrust  response  is  acceptable,  EGT  appears  to  be  uncontrolled.  The  need  for  gain  scheduling  of  robust 
controllers  is  clearly  evident  due  to  the  nonlinear  effects  of  engine  deterioration. 


4 


* 
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Figure  5.5. — 70%  degraded  engine  response  with  LPV  nominal  engine  controller 
a)  Thrust  response,  b)  EGT  deviation  response,  c)  VSV  and  VBV  response. 
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An  example  of  controller  robustness  to  engine-to-engine  variation  is  shown  in  Figure  5.6.  The  base  values  of  the 
health  parameters  correspond  to  70%  engine  deterioration  -  which  are  taken  to  be  fleet  averages.  A  Monte  Carlo 
approach,  using  standard  deviations  calculated  from  fleet  data,  is  used  to  determine  a  health  parameter  set  for  each 
specific  engine.  As  with  the  example  presented  in  Figures  5.1  and  5.2  (the  40%  deteriorated  case),  the  net  effect  is  a 
steady-state  reduction  in  EGT.  For  the  engines  considered,  EGT  steady-state  reductions  range  from  2  to  4  degrees, 
with  the  corresponding  thrust  response  virtually  identical  for  each  engine.  Compressor  stall  margin  deviation  is  also 
shown  to  illustrate  that  component  margins  remain  well  above  their  critical  values  for  the  engines  considered. 


Time  dec) 

Figure  5.6.— Deteriorated  engine  response  with  LPV  controller  showing 
robustness  to  engine-to-engine  variation  about  an  EDI  of  0.7  (dashed  line), 
a)  Thrust  response,  b)  compressor  SMD  response,  c)  EGT  deviation  response. 


VI.  Conclusions 

An  initial  study,  via  testing  with  a  nonlinear  engine  simulation,  has  demonstrated  the  ability  of  LPVQLF-based 
control  to  simultaneously  control  engine  thrust  and  exhaust  gas  temperature  while  being  robust  to  engine 
deterioration  and  engine-to-engine  variation.  By  reducing  EGT  for  a  given  thrust  level,  it  might  be  possible  to 
reduce  TMF  damage  and  thus  extend  the  on-wing  life  of  the  engine.  An  LMI-based  optimization,  searching  for 
feasible  solutions  to  quadratic  Lyapunov  functions  based  on  engine  state  space  dynamics,  provides  a  means  for 
designing  unconditionally  stable  gain  scheduled  controllers  across  a  full  range  of  deteriorated  engines.  As  well,  the 
quadratic  Lyapunov  function  approach  guarantees  stability  for  instantaneous  scheduling  parameter  variations. 
Traditionally  gain  scheduling  is  based  on  an  ad  hoc  approach,  requiring  multiple  pre-designed  LTI  controllers. 
These  techniques,  while  employing  controllers  shown  to  be  stable  at  the  points  of  plant  linearization,  are  not 
guaranteed  to  be  stable  for  plants  between  those  for  which  the  controllers  were  specifically  designed. 

With  the  LPVQLF  approach,  system  (i.e.,  engine-to-engine)  uncertainty  is  handled  by  operating  condition 
dependant  weights  while  other,  more  significant  uncertainties  due  to  the  highly  nonlinear  nature  of  the  plant,  are 
addressed  by  use  of  the  gain  scheduled  LPV  controller.  The  control  system  is  designed,  for  cruise  conditions,  using  a 
feedforward  neural  network  to  reduce  the  set  of  estimated  health  parameters  to  a  single  parameter  for  gain 
scheduling,  the  engine  deterioration  index  EDI.  Future  work  may  include  controller  design  over  a  larger  region  of 
the  flight  envelope,  possibly  targeting  more  critical  phases  of  operation  such  as  the  climb  phase  of  flight.  In  addition, 
as  is  the  case  in  many  robust  controller  designs,3'6  the  models  used  to  synthesize  the  controllers  were  not  pre-scaled. 
Future  work  may  also  investigate  pre-scaling  as  a  means  to  enhance  the  performance  of  the  LPV  controller. 
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The  LPV  controller  is  shown  to  be  stable  for  all  scheduling  parameter  trajectories,  but  the  trajectories  are  still  based 
on  a  nominal  degradation  profile.  Although  the  controller  is  robust  to  wide  health  parameter  variation  around  the 
nominal  trajectory,  it  does  not  take  into  account  large  off-nominal  deviations  of  individual  health  parameters. 
Likewise,  the  health-parameter-to-EDI  mapping  employed  (i.e.,  a  feedforward  neural  network)  was  developed  using 
the  nominal  degradation  profile.  Values  of  EDI  produced  by  the  neural  network  as  a  result  of  off-nominal  health 
parameter  input  will  be  invalid,  potentially  leading  to  the  scheduling  of  an  inappropriate  controller.  Thus  future 
work  should  look  at  the  issue  robustness  of  the  EDI  mapping  to  off-nominal  health  parameter  variation,  and  its 
relationship  to  the  scheduled  controller  to  understand  the  limitations  imposed  by  the  EDI  on  the  effectiveness  of  the 
LPV  controller  implementation. 

Due  to  the  use  of  quadratic  Lyapunov  functions  and  the  corresponding  guaranteed  stability  for  instantaneous 
scheduling  parameter  variations,  combining  the  technique  described  with  a  method  for  assessing  engine  condition 
e.g.,  use  of  finite  state  machines  or  discrete  event  systems13  to  provide  a  switching  signal,  may  prove  to  be  extremely 
powerful  for  implementing  supervisory  robust  control  of  engines  that  incur  damage  -  especially  for  those  damage 
scenarios  that  may  not  be  fully  (or  correctly)  hypothesized. 
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